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Classical rotational inertia of solid 4 He 
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The observation of reduced rotational inertia in a cell containing solid 4 He has been interpreted 
as evidence for superfluidity of the solid. An alternative explanation is slippage of the solid at the 
container wall due to grain boundary premelting between the solid and dense adsorbed layers at 
the container wall. We calculate the range of film thickness and the viscous drag, and find that the 
slippage can account for the observations. 

PACS numbers: 67.40.Hf, 67.80.-s, 67.70.+n, 67.57.De, 67.57.Np 



A recent report by Kim and Chan [l| describes the 
observation of inertial anomalies, termed "non classical 
rotational inertial fraction (NCRIF)" in solid 4 He, which 
are taken to demonstrate superfluidity of the solid. The 
exciting possibilities raised by the experiment and their 
earlier study Q wherein 4 He was solidified in Vycor, a 
microporous glass, span a range of fundamental issues in 
quantum materials. Leggett notes the possibility of 
supersolid behavior as a long standing speculation, based 
on the hypothesis of Bose-Einstein condensation of zero- 
point vacancies, and Prokof'ev and Svistunov argue |4| 
that the vacancy density in well ordered solid He is insuf- 
ficient for Bose-Einstein condensation at the experimen- 
tal temperatures, speculating instead that the NCRIF in 
the Vycor study may be due to the superfluidity of 
vacancies in a defect-laden layer of solid 4 He at the Vy- 
cor surfaces. Finally, Beamish [5j noted the possibility 
that a disordered thin liquid like layer at the Vycor walls 
may persist at low temperature and argued that its su- 
perfluid properties would be different from those found 
at pressures well below solidification pressure. 

Here we suggest an alternative explanation of the ex- 
periment on bulk solid [1]: slippage of the solid at the 
wall of the container due to a liquid film caused by grain 
boundary premelting. The premelting in this case is not 
at ordinary grain boundaries, but at the interface be- 
tween the bulk solid and dense adsorbed layers at the 
container wall. The dense layers, due to strong adsorp- 
tion forces, are responsible for nonzero wetting angles be- 
tween solid 4 He and copper and glass walls 0, Q ; rather 
than an ordinary grain boundary, the contacting surface 
in question more nearly resembles the interface between 
two different materials. The liquid film separates the 
bulk solid from the torsional balance and replaces the 
shear strength of the dry interface by viscous drag. Our 
model allows the possibility that the film is superfluid, in 
which case the calculation is equally applicable to its nor- 
mal fluid component. In the following we first estimate 
the thickness of the premelted film, and then calculate 
the dynamical coupling. 

The nature of premelting in any system is determined 
by the competition between bulk and surface free ener- 
gies. Complete premelting, in which the thickness of the 
melt layer diverges as temperature approaches the bulk 



transition, requires that the total excess surface free en- 
ergy per unit area, F(d), be a positive monotonically 
decreasing function of the film thickness with a global 
minimum at infinite film thickness. In incomplete pre- 
melting, the melt thickness remains finite at the bulk 
transition, such as in recent studies of ice It is a 
general result 0-01 that in a symmetrical system (e.g., 
solid/liquid/solid) the long range interactions are attrac- 
tive, and consequently grain boundary melting must be 
incomplete. The most complete formulation of the excess 
surface free energy per unit area, F(d), for systems en- 
tirely controlled by frequency dependent dispersion forces 
is that of Dyzaloshinskii, Lifshitz and Pitaevskii (DLP) 
|12| . However, it requires as input the frequency de- 
pendent dielectric properties of the layers in the system 
under consideration. Because we are dealing with bulk 
solid 4 He and the dense adsorbed solid at the container 
wall, the input data for the DLP theory are not avail- 
able. Therefore, we proceed with reasonable ranges of 
the Lennard- Jones parameters 0, Q] . 

The total free energy of the system at a given temper- 
ature T and pressure P is written 

G T (T,P,d) = Pm (T,P)d + F(d), (1) 

where the liquid density and chemical potential are pi 
and fii(T, P). In grain boundary premelting 
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F{d)= lss {d) = 2 lsl + Pl / V(z)dz, (2) 



where 7 ss (ci) is the interfacial free energy of the solid- 
solid interface, and j s i is the interfacial free energy per 
unit area of the solid-liquid interfaces, with implicit ref- 
erence to the crystallographic orientation present at an 
interface. In lieu of the DLP theory, the most general 
phenomenological mean field model considers V(z) as the 
Lennard-Jones potential 13] but augmented to embody 
the effects of retardation viz a viz 
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where z is the coordinate normal to the surfaces, C3 (B) 
is the nonretarded (retarded) van der Waals attraction, 
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and D is the well depth. At each temperature and pres- 
sure below T\ the bulk and intcrfacial free energies strike 
a balance and one can show [e.g., 0,0] that a unique 
equilibrium film thickness obtains from 
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where the latent heat of fusion is q m and t is the re- 
duced temperature. From this we can simply write 
the equilibrium film thickness-temperature relation as 
t = q m - X V{d). 

In the figure we plot d = d(t) for a range of potential 
parameters suggested from the detailed analysis of the 
wetting of wide classes of substrates by liquid helium 0] . 
Curves for three values of the latent heat of fusion q m are 
shown and the value at the NCRIF onset temperature of 
175 mK is extrapolated to lower temperatures from 1 K 
using the data of Swenson [lfi| . Because t = q m V(d), 
we find q m to be the most important parameter in the 
problem: The latent heat embodies the bulk free energy 
penalty for converting solid to liquid, against the melt 
driving interactions of the potential. Therefore, the gen- 
eral understanding that q m — > as T — > (t — > 1), but 
the lack of experimental information on q rn , at solidifica- 
tion pressures and temperatures below about IK, leaves 
open important quantitative questions. Thus, we view 
our film thickness calculations as conservative-thinner 
than is likely. We find that although the magnitude of 
the thickness of a premelted layer depends on the param- 
eters used in the calculation, (a) the temperature depen- 
dence itself is rather weak and (b) the film thicknesses 
(1-4 atomic layers) are sufficient to accommodate flow 
and superflow (e.g., |l7LIT^). 

We now estimate the dynamical coupling of the solid 
to the container, through the premelted layer. The solid 
is a thin walled hollow cylinder, of mean radius a, wall 
thickness s and height h, bathed on both sides by a 
layer of liquid thickness d. The liquid is driven by the 
container's torsional oscillations of angular displacement 



6{t) = 



driving the solid at the same frequency, 



and angular amplitude 0' o . The hydrodynamic regime is 
governed by the relative magnitudes of d and the decay 
length of transverse viscous waves A = y/2rj/piU>, where 
the dynamic viscosity is rj and the density of the liquid 
19, 201. To evaluate A we have the experimental fre- 



2.3 x 10 3 s -1 and the liquid density (l(i| 
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pi w 0.2g/cm in the range of experimental pressures. 
Since the viscosity of the premelted liquid is not known, 
we must estimate it from the properties of the bulk liq- 
uid. The viscosity of normal liquid Hel is about 20p,P; 
the viscosity of the normal component of Hell between 
1.2 K and T\ ranges between 1(VP and 20^P 20]. How- 
ever, the viscosity may be enhanced in the very narrow 
gap d [T^i thereby placing an uncertainty on the exper- 
imental value of rj. As a conservative measure, we will 
allow a possible enhancement of an order of magnitude. 
Accordingly, we estimate 2 x 10 -4 < A < 7 x 10~ 4 cm, 
which is on the order of one thousand times the thickness 



d of the grain boundary melted liquid, thereby reducing 
the problem to one of slow, nearly steady flow. A second 
simplification stems from the ratio of d to the cylinder 
radius a; d/a << 1, which makes the problem equivalent 
to the drag between parallel plates. Therefore the fluid 
velocity varies linearly between the surface of the cell and 
the solid helium, so that the viscous drag per unit area 
on the solid helium is / = (aojrj/d)(8 - 6' )e iuJt . The 
torque Tt due to the total force on both inner and outer 
surfaces of the solid cylinder is 



4ira 3 hujri 
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The torque induces the solid's inertial response, the time 
rate of change of angular momentum: L s = Iu 2 0' o e wt , 
where I — 2ira 3 hsp so i, in which we emphasize that p so i is 
the density of the solid to avoid confusion with the con- 
ventional nomenclature for the superfluid density. Thus, 
the fractional difference in the amplitudes of rotational 
motion between the bulk solid and the premelted film is 
written as 6 ° g , ° = S, and hence the fractional inertial re- 
sponse is e ° e 9 ° — , wherein the controlling parameter 
is the value of S which is written as follows 

j ds Psol 

A^ pi 



(6) 



The ranges of possible values of d and r\ impose wide 
limits on the estimate of slippage: 0.004 < 6 < 0.26. 
These extremes span the maximum values of the NCRIF 
measured at low temperature and low amplitude, roughly 
0.005 to 0.02, depending on the pressure 0. 

Additional observations of Kim and Chan are that 
the NCRIF decreases at higher amplitudes of oscillation, 
which is given as strong supporting evidence of superfluid 
solid behavior, attributing it to exceeding the critical ve- 
locity of superfluidity in the solid . We propose that it 
can indeed be evidence of superfluidity, but in the liquid 
film rather than in the solid. The slippage is controlled, 
at low relative speeds, by the viscosity of the normal 
fluid component, and is augmented, at higher speed, by 
the excitation of the superfluid fraction. 

The results of the earlier work with Vycor are rele- 
vant here. Because of pinning of the solid by the Vycor 
matrix Q , the loss of rotational inertia is interpreted as 
an indication of superfluidity It may be, as Prokof 'ev 
and Svistunov 4] and Beamish suggest, in a disor- 
dered solid layer at the Vycor interface. But it can be 
expected that grain boundary melting occurs in the Vy- 
cor study, as in the experiment with bulk solid. Thus 
the superfluidity would be associated with the liquid film, 
rather than the solid. Indeed, the nm scale disorder of 
Vycor was used in order to enhance the density of va- 
cancies 0), and depending on pressure and the assumed 
fraction of material locked by the tortuosity of the Vycor, 
the estimate of the BEC fraction could be less than that 
in the bulk solid. In the framework of our theory, the 
difference between the experiments may be explained in 
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terms of the wall materials inducing a different density of 
adsorbed solid. Interesting tests that span the existing 
range of behavior would involve the insertion of a ridge in 
the cell wall of the most recent experiment or changing 
the wall material. An appropriate ridge dimension would 
lock the solid and thereby offer a possible test of whether 
slippage had explained the NCRIF. An additional obser- 
vation is an extreme sensitivity of the NCRIF to 3 He; in 
our model it is explained by the very small fraction of 
the sample that is involved in the grain boundary melt, 



and the greater solubility of 3 He in liquid, rather than in 
solid helium. 

Finally, although we propose an alternative explana- 
tion to a superfiuid solid, we consider that decoupling 
due to a premelted, and possibly superfiuid, film offers 
interesting possibilities. The film exists in a region of 
pressure and temperature not otherwise accessible, and 
the experiment suggests interesting and entirely new pos- 
sibilities in studies of premelting and liquid confinement. 
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FIG. 1: Plots of liquid 4 He film thickness, in units of atomic 
layers, as a function of reduced temperature, t, and values of 
the Lennard- Jones potential parameters C3/D = 102.9, 55.3, 
and 12.5 A (as labeled adjacent to the curves) chosen for 
expected ranges fl3T | and because retardation typically occurs 
at ranges that are much larger than is found, we let B =0. The 
two different sets of three curves are meant to demonstrate 
the important influence of the latent heat of fusion, q m ; (a) 
q m (T = 1.2K) = 0.09 cal/mole, and (b) q m (T = 175mK) = 
0.0004 cal/mole. The data from Swenson [lg] are extrapolated 
to the experimental temperature of 175 mK based on the 
requirement that q m vanishes at absolute zero. 



